Large deformed structures in A^e — S nuclei near neutron drip-line 
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The structure of Ne, Na, Mg, Al, Si and S nuclei near the neutron drip-line region is investigated in the frame- 
work of relativistic meannfield (RMF) and non-relativistic Skyrme Hartree-Fock formalisms. The drip-line of 
these nuclei are pointed out. We analysed the large deformation structures and many of these neutron rich nuclei 
are quite deformed. New magic number are seen for these deformed nuclei. 
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I. INTRODUCTION 

The structure of light nuclei near the neutron drip-line is 
very interesting for a good number of exotic phenomena. Nu- 
clei in this region are very different in collectivity and clus- 
tering features than the stable counterpart in the nuclear chart. 
For example, the neutron magicity is lost for the N=8 nucleus 
for i^Be [1] and N=20 for ^^Mg [2]. The discovery of large 
collectivity of '^"'Mg by Iwasaki et al. [3] is another example of 
such properties. The deformed structures and core excitations 
of Mg and neighboring nuclei and location of drip-line in this 
mass region is an important matter 0]. On the other hand, the 
appearance of N=16 magic number for ^'^O is well established 
15]. The discovery of the two isotopes "'"Mg and "^^Al, once 
predicted to be drip-line nuclei is |7|] gives indication that the 
neutron drip-line is located towards the heavier mass region. 
The existence of neutron halo in ^^Li is well established and 
the possibility of proton halo in '^B and the neutron halo in 
^^Be and ^^B are very interesting phenomena for the drip-line 
nuclei. In addition to these above exciting properties, the clus- 
ter structure of entire light mass nuclei and the skin formation 
in neutron-drip nuclei provide us features for the study of light 
mass drip-line nuclei. Also, the exotic neutron drip-line nuclei 
play a role in many astrophysical studies. In this paper, our 
aim is to study the neutron drip-line for Ne-S isotopic chain 
in the frame-work of a relativistic mean field (RMF) and non- 
relativistic Skyrme Hartree-Fock formalism and analyse the 
large deformation of these isotopes. 

The paper is organised as follows: The relativistic and non- 
relativistic mean field formalisms are described very briefly 
in Section II. The results obtained from the relativistic mean 
field (RMF) and Skyrme-Hartree-Fock(SHF) formalisms, and 
a discussion of these results, are presented in Section III. Fi- 
nally summary and concluding remarks are given in Section 
IV. 



II. THEORETICAL FRAMEWORK 

Mean field methods have been widely used in the study of 
binding energies and other properties of nuclei ^ . Al- 
though the older version of the SHF and RMF models have 
some limitation to reproduce some of the observables, the re- 
cent formalisms are quite efficient to predict the bulk prop- 
erties of nuclei not only near the stability valley, but also for 



the nuclei near the proton and neutron drip-Unes. We use here 
two of the successful mean field models [Tol - [T6ll (Skyrme 
Hartree-Fock and the Relativistic Mean Field ) to learn about 
the properties of drip-line nuclei A^e — S. 



A. The Skyrme Hartree-Fock (SHF) Method 

There are many known parametrizations of Skyrme interac- 
tion which reproduce the experimental data for ground-state 
properties of finite nuclei and for the observables of infinite 
nuclear matter at saturation densities, giving more or less 
comparable agreements with the experimental or expected 
empirical data. The general form of the Skyrme effective in- 
teraction, used in the mean-field models, can be expressed as 
a density functional liKXill.l . given as a function of some 
empirical parameters, as 



n = ic + no + ns + n 



eff 



(1) 



where JC is the kinetic energy term, Ho the zero range, TLs 
the density dependent and "He// the effective-mass depen- 
dent terms, which are relevant for calculating the properties 
of nuclear matter. These are functions of 9 parameters ti, Xi 
{i — 0,1,2,3) and 77, and are given as 



-to [(2 + xo)p^ - {2x0 + 1){pI + pD] , (2) 



Tio 

^3 = :^hp^[{2 + x3)p' ~{2x3 + 1){pI + pI)], (3) 
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= l[ti{2 + Xi)+t2{2+X2)]Tp 



+ - [t2i2x2 + 1) - h{2xi + 1)] {Tppp + TnPn)- 
o 

(4) 



The kinetic energy K, ~ ^'t^ ^ form used in the Fermi gas 
model for non-interacting fermions. The other terms, repre- 
senting the surface contributions of a finite nucleus with 64 



2 



and 64 as additional parameters, are 
1 
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1, 



and 



bipV ■ J + b'^{pnV ■ Jn + PpV ■ Jp) 



(5) 
.(6) 



Here, the total nucleon number density p — Pn + Pp, the ki- 
netic energy density r = t„ + r^, and the spin-orbit density 
J = Jn + Jp- The subscripts n and p refer to neutron and 
proton, respectively, and rn is the nucleon mass. The Jq — 0, 
q = n OT p, for spin-saturated nuclei, i.e., for nuclei with ma- 
jor oscillator shells completely filled. The total binding energy 
(BE) of a nucleus is the integral of the density functional H. 

At least eighty-seven parametrizations of the Skyrme inter- 
action are published since 1972 lfT2ll where 64 = 64 = Wq, we 
have used here the Skyrme SkI4 set with 64 ^ 64 fo*]. This 
parameter set is designed for considerations of proper spin- 
orbit interaction in finite nuclei, related to the isotope shifts 
in Pb region and is better suited for the study of exotic nu- 
clei. Several more recent Skyrme parameters such as SLyl- 
10, SkX, SkI5 and SkI6 are obtained by fitting the Hartree- 
Fock (HF) results with experimental data for nuclei starting 
from the valley of stability to neutron and proton drip-lines 



rrom the valley 



I (41^~^/^). The constant gap BCS pairing is used to add the 
pairing effects for the open shell nuclei. It is to be noted that 
in the present work only intrinsic state solutions are presented. 
Each of these deformed intrinsic states is a superposition of 
various angular momenta states. To obtain the good angu- 
lar momentum states and spectroscopic predictions for these 
nuclei near neutron drip-line we need to project out states of 
good angular momenta. Such calculation will be considered 
as a future extension of this work. The quadrupole moment 
deformation parameter (32 is evaluated from the resulting pro- 
ton and neutron quadrupole moments, as Q = Qn + Qp = 

\J ^^{-^AR^I32)- The root mean square (rms) matter ra- 
dius is defined as (r^j) = ^/ p{r±, z)r'^dT; here A is the 
mass number, and p{r±, z) is the deformed density. The to- 
tal binding energy and other observables are also obtained by 
using the standard relations, given in iietl . We use the well 
known NL3 parameter set |i 1 8i1 . This set not only reproduces 
the properties of stable nuclei but also well predicts for those 
far from the /^-stability valley. Also, the isoscalar monopole 
energy agrees excellently with the experimental values for dif- 
ferent regions of the Periodic Table. The measured superde- 
formed minimum in ^^''Hg is 6.02 MeV above the ground 
state, whereas in RMF calculation with NL3 set, this num- 
ber is 5.99 MeV |18]. All these facts give us confidence to 
use this older, though very much still in use, NL3 set for the 
present investigation. 



m. RESULTS AND DISCUSSION 



B. The Relativistic Mean Field (RMF) Method 

The relativistic mean field approach is well-known and the 
theory is well documented lll5lll6ll . Here we start with the rel- 
ativistic Lagrangian density for a nucleon-meson many-body 
system, as 

C = W^{ird^, - M}V, + ^d^'ad^a - 

-If^^F,,, - eV^,7^^^^V'.^M• (7) 

All the quantities have their usual well known meanings. 
From the relativistic Lagrangian we obtain the field equations 
for the nucleons and mesons. These equations are solved 
by expanding the upper and lower components of the Dirac 
spinors and the boson fields in an axially deformed harmonic 
oscillator basis with an initial deformation. The set of coupled 
equations is solved numerically by a self-consistent iteration 
method. 

The centre-of-mass motion energy correction is esti- 
mated by the usual harmonic oscillator formula Ec.m. = 



A. Ground state properties from the SHF and RMF models 

There exists a number of parameter sets for the standard 
SHF and RMF Hamiltonians and Lagrangians. In some of our 
previous papers and of other authors I Iq, [T8I - I21I1 the ground 
state properties, like the binding energies (BE), quadrupole 
moment deformation parameters /32, charge radii {tc) and 
other bulk properties are evaluated by using the various non- 
relativistic and relativistic parameter sets. It is found that, 
more or less, most of the recent parameters reproduce well 
the ground state properties, not only of stable normal nuclei 
but also of exotic nuclei which are away from the valley of 
beta-stability. So, if one uses a reasonably acceptable param- 
eter set the prediction of the results will remain nearly force 
independent. This is valid both for SHF and RMF formalisms. 
However, with a careful inspection of these parametrizations, 
some of the SHF and RMF sets can not reproduce the empir- 
ical data. In this context we can cite the deviation of isotopic 
shifts than the experimental data |22] for Pb nuclei while us- 
ing SHF forces like, SkM* values 112311 . However, the RMF 
sets reproduce the kink quite nicely On the other hand, 
most of the RMF sets over estimate the nuclear matter incom- 
pressibility. In general, the predictive power of both the for- 
malisms are reasonably well and can be comparable to each 
other, which can be seen in the subsequent subsections. In ad- 
dition to this, the general results SHF (SkM) and RMF (NL3) 
forces are similar for the considered region. Thus in the sub- 
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TABLE I: The calculated ground state binding energy obtained from 
SHF and RMF theory are compared with the experimentally known 
heaviest isotope for Ne, Na, Mg, Al, Si and S n23]. 



nucleus 


RMF 


SHF 


Expt. 


nucleus 


RMF 


SHF 


Expt. 




215.1 
257.7 
310.1 


210.6 
255.1 
305.1 


211.2 
256.2 
306.6 


33 Na 

39^1 
45 s 


237.9 
285.2 
353.5 


234.5 
281.9 
350.2 


232.8 
283.1 
354.2 



TABLE II: The predicted neutron drip-line for Ne, Na, Mg, Al, Si 
and S in RMF (NL3) and SHF (SKI4) parameter sets are compared 
with prediction of infinite nuclear matter (INM) mass model ^2^1 , 
finite range droplet model (FRDM) 1,27.] and experimental data i23l/ 
along with the number shown in parenthesis are the experimentally 
extrapolated values. 



nucleus 


RMF 


SHF 


INM 


FRDM 


Expt. 


Ne 


34 


34 


36 


32 


30 (34) 


Na 


40 


37 


40 


36 


33 (37) 


Mg 


40 


40 


46 


40 


34 (40) 


Al 


48 


48 


48 


45 


39 (42) 


Si 


54 


48 


50 


48 


41 (44) 


S 


55 


55 


53 


50 


45 (49) 



sequent results during our discussion we will refer the results 
of RMF (NL3) calculations, except some specific cases. Thus 
the result of SHF (SkI4) are not displayed in Tables. 



B. Binding energy and neutron drip-line 

The ground state binding energy (BE) are calculated for Ne, 
Na, Mg, Al, Si and S isotopes near the neutron drip-line. This 
is done by comparing the prolate, oblate and spherical solu- 
tion of binding energy for a particular nucleus. For a given nu- 
cleus, the maximum binding energy corresponds to the ground 
state and other solutions are obtained as various excited intrin- 
sic states. In Table I, the ground state binding energy for the 
heaviest isotopes for the nuclei discussed are compared with 
the experimental data |25]. From the Table it is observed that 
the calculated binding energies are comparable with SHF and 
RMF results. We have listed the neutron drip-lines in Table 
II, which are obtained from the ground state binding energy 
for neutron rich Ne, Na, Mg, Al, Si and S nuclei. The nu- 
clei with the largest neutron numbers so far experimentally 
detected in an isotopic chain till date, known as experimental 
neutron drip-line are also displayed in this Table for compari- 
son. The numbers given in the parenthesis are the experimen- 
tally extrapolated values J25t]. To get a qualitative understand- 
ing of the prediction of neutron drip-line, we have compared 
our results with the infinite nuclear matter (INM) [26.1 and fi- 
nite range droplet model (FRDM) ll27ll mass estimation. From 
the table, it is clear that all the predictions for neutron drip-line 
are comparable to each other. 

The drip-lines are very important after discovery of the two 
isotopes '"^Mg and '^^Al |6|1 that here once predicted to be be- 
yond the drip-line 10, l28ll . This suggests that the drip-line is 
somewhere in the heavier side of the mass prediction which 
are beyond the scope of the present mass models ll7l l28ll . In 



TABLE III: The calculated value of charge radius (Tch), quadrupole 
moment deformation parameter /32 arid binding energy (BE) for Ne, 
Na and Mg nuclei in RMF (NL3) formalism. The maximum binding 
energy is the ground state solution and all other values are the in- 
trinsic excited state solution. The radius Teh is infm and the binding 
energy is in MeV. 



nucleus 


Teh 


132 


BE(MeV) 


nucleus 


rch 


P2 


BE(MeV) 


""Ne 


2.970 


0.535 


156.7 


""Ne 


2.901 


-0.244 


152.0 


2iNe 


2.953 


0.516 


165.9 


^^Ne 


2.889 


-0.241 


161.1 


^^Ne 


2.940 


0.502 


175.7 


^^Ne 


2.881 


-0.242 


170.5 


^^'Ne 


2.913 


0.386 


181.8 


2^Ne 


2.880 


-0.249 


179.7 


^''Ne 


2.890 


0.278 


188.9 


2*Ne 


2.881 


-0.259 


188.9 


^^Ne 


2.907 


0.272 


194.2 


^^Ne 


2.886 


-0.206 


194.2 


2«Ne 


2.926 


0.277 


199.9 


26Ne 


2.893 


-0.159 


199.9 


^■'Ne 


2.945 


0.247 


203.9 


^''Ne 


2.925 


-0.183 


203.9 


2»Ne 


2.965 


0.225 


208.2 


28Ne 


2.957 


-0.203 


208.2 


29Ne 


2.981 


0.161 


211.2 


29Ne 


2.974 


-0.133 


211.2 


3°Ne 


2.998 


0.100 


215.0 


^''Ne 


2.995 


-0.081 


215.0 


3^Ne 


3.031 


0.244 


216.0 


^^Ne 


3.013 


-0.133 


216.0 


^^Ne 


3.071 


0.373 


218.6 


^^Ne 


3.033 


-0.180 


218.6 


^^'Ne 


3.095 


0.424 


219.5 


3^Ne 


3.044 


-0.230 


219.5 


34Ne 


3.119 


0.473 


220.9 


^^Ne 


3.054 


-0.275 


220.9 


^'^Ne 


3.132 


0.505 


220.4 


^^Ne 


3.064 


-0.315 


215.7 


3«Ne 


3.146 


0.539 


220.3 


3'^Ne 


3.075 


-0.352 


220.4 




2.939 


- 0.250 


189.4 


^''Na 


2.964 


0.379 


192.3 


^^Na 


2.938 


-0.258 


200.3 


25Na 


2.937 


0.273 


200.6 


2«Na 


2.940 


-0.202 


206.3 


2«Na 


2.965 


0.295 


207.1 


^■^Na 


2.946 


-0.157 


212.5 


^^Na 


2.993 


0.323 


214.2 


2«Na 


2.980 


-0.184 


217.7 


2«Na 


2.993 


0.272 


219.0 


^^Na 


3.012 


-0.205 


223.4 


^'■'Na 


3.004 


0.232 


224.3 


3°Na 


3.025 


-0.131 


227.5 


^''Na 


3.031 


0.169 


228.1 


^^Na 


3.043 


-0.074 


232.5 


^^Na 


3.047 


0.108 


232.7 


^^Na 


3.061 


-0.129 


233.0 


^^Na 


3.077 


0.237 


234.5 


^^Na 


3.082 


-0.179 


234.3 


^^Na 


3.113 


0.356 


237.9 


^''Na 


3.095 


-0.226 


234.8 


3*Na 


3.137 


0.404 


239.8 


^^Na 


3.108 


-0.270 


235.9 


^^Na 


3.161 


0.450 


242.3 


^^Na 


3.121 


-0.308 


236.9 


3'^Na 


3.175 


0.481 


242.5 


^■^Na 


3.135 


-0.345 


238.4 


3^Na 


3.190 


0.512 


243.1 


3«Na 


3.156 


-0.359 


240.0 


38Na 


3.199 


0.491 


243.4 


3^Na 


3.180 


-0.375 


241.8 


^''Na 


3.209 


0.472 


244.1 


«Na 


3.184 


-0.358 


241.3 


*''Na 


3.228 


0.477 


243.4 


'"Mg 


3.043 


0.487 


194.3 


'"Mg 


3.001 


-0.256 


186.8 


2^Mg 


3.009 


0.376 


202.9 


^^Mg 


2.993 


-0.261 


199.1 


^''Mg 


2.978 


0.273 


212.5 


^^Mg 


2.990 


-0.268 


211.6 


''Mg 


3.015 


0.310 


220.2 


''Mg 


2.988 


-0.204 


218.2 


2**Mg 


3.048 


0.345 


228.7 


2»Mg 


2.992 


-0.154 


225.6 


2«Mg 


3.055 


0.289 


234.3 


29Mg 


3.027 


-0.186 


232.0 




3.062 


0.241 


240.5 


»°Mg 


3.059 


-0.207 


239.0 


30 Mg 


3.131 


0.599 


237.7 










^^Mg 


3.075 


0.179 


245.1 


3iMg 


3.068 


-0.128 


244.0 


^'^Mg 


3.090 


0.119 


250.5 


^^Mg 


3.085 


-0.067 


249.9 


^^Mg 


3.131 


0.471 


248.8 










33 Mg 


3.117 


0.233 


253.1 


33Mg 


3.102 


-0.126 


251.5 


^■^Mg 


3.150 


0.343 


257.3 


»4Mg 


3.124 


-0.181 


253.9 


^*Mg 


3.184 


0.588 


254.1 










35 Mg 


3.173 


0.388 


260.5 


35Mg 


3.141 


-0.196 


255.3 


36 Mg 


3.198 


0.432 


263.9 


36Mg 


3.160 


-0.213 


257.0 


^^Mg 


3.212 


0.462 


264.9 


^■^Mg 


3.179 


-0.258 


258.8 


38 Mg 


3.227 


0.492 


266.3 


38Mg 


3.198 


-0.300 


261.1 


39 Mg 


3.237 


0.473 


267.8 


39Mg 


3.216 


-0.338 


263.4 


«Mg 


3.247 


0.456 


269.7 


«Mg 


3.234 


-0.374 


266.4 
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TABLE IV: Same as Table III, but for Al, Si and S. 



nucleus 


Teh 


/32 


BE(MeV) 


nucleus 


Teh 


02 


BE(MeV) 




3.097 


0.388 


182.3 


^''Al 


3.077 


-0.258 


179.4 


^^Al 


3.072 


0.381 


197.7 


^-''Al 


3.060 


-0.266 


193.9 


26 


3.191 


0.550 


206.6 


26 


3.052 


-0.275 


207.8 


^■^Al 


3.215 


0.572 


217.0 


2^A1 


3.053 


-0.292 


221.9 


2*A1 


3.178 


0.471 


226.7 


28 


3.037 


-0.208 


238.6 


^®A1 


3.061 


0.251 


239.3 


2"A1 


3.033 


-0.141 


245.6 


^°A1 


3.073 


0.207 


246.2 


3"A1 


3.070 


-0.184 


253.8 


31 Al 


3.085 


0.170 


253.6 


31 Al 


3.101 


-0.205 


259.8 


32 Al 


3.101 


0.113 


260.0 


32 Al 


3.103 


-0.111 


261.2 


33 Al 


3.118 


0.057 


267.2 


33 Al 


3.165 


-0.333 


269.4 


3''A1 


3.139 


0.159 


269.9 


3"A1 


3.134 


-0.108 


275.1 


35 Al 


3.167 


0.268 


274.1 


35 


3.157 


-0.172 


272.8 




3.187 


0.313 


277.6 


3«A1 


3.173 


-0.189 


277.7 


3'^Al 


3.208 


0.355 


281.5 


3^A1 


3.191 


-0.208 


280.3 


3»A1 


3.275 


0.418 


282.7 


3^A1 


3.214 


-0.254 


283.5 


3^A1 


3.285 


0.406 


285.1 


3'^Al 


3.236 


-0.299 


286.7 


^°A1 


3.304 


0.441 


287.6 


""Al 


3.257 


-0.336 


290.4 


''^Al 


3.325 


0.474 


290.5 


"^Al 


3.278 


-0.370 


290.6 


"'^Al 


3.348 


0.483 


290.7 


"^Al 


3.281 


-0.355 


291.2 


''3A1 


3.371 


0.491 


291.3 


"3A1 


3.282 


-0.338 


292.2 


""Al 


3.375 


0.456 


291.0 


""Al 


3.274 


-0.288 


293.6 


"^Al 


3.378 


0.420 


291.0 


"■■"'Al 


3.271 


-0.263 


293.5 


"•'Al 


3.359 


0.341 


294.5 


"*'A1 


3.346 


-0.296 


294.0 


'"'Al 


3.246 


0.125 


293.6 


"*'A1 


3.432 


0.660 


290.8 


"■^Al 


3.246 


0.090 


294.8 


"^Al 


3.335 


-0.319 


293.6 


"■^Al 


3.447 


0.653 


290.4 










"*A1 


3.276 


0.117 


293.8 


"**A1 


3.319 


-0.252 


294.0 


2" Si 


3.162 


0.294 


169.1 


2"Si 


3.170 


-0.278 


169.3 


2^81 


3.127 


0.286 


185.1 


25 Si 


3.139 


-0.274 


184.8 




3.099 


0.282 


201.8 


Si 


3.118 


-0.280 


200.9 


2^Si 


3.054 


0.168 


215.8 


2''Si 


3.114 


-0.299 


216.4 


^®Si 


3.017 


0.001 


231.4 


2** Si 


3.122 


-0.331 


232.1 




3.035 


0.001 


240.7 


2''Si 


3.093 


-0.237 


240.7 


30 Si 


3.070 


0.148 


250.6 


3° Si 


3.054 


-0.060 


250.4 


3iSi 


3.089 


0.120 


258.7 


31 Si 


3.108 


-0.180 


259.1 


32si 


3.109 


0.104 


267.2 


32 Si 


3.137 


-0.201 


268.5 


33 Si 


3.126 


0.050 


275.4 


33 Si 


3.131 


-0.084 


275.6 


3" Si 


3.148 


0.000 


284.4 


3" Si 


3.204 


-0.336 


278.5 


35 Si 


3.160 


0.085 


287.3 


35 Si 


3.161 


-0.083 


287.4 


3«'Si 


3.184 


0.193 


291.4 


36 Si 


3.186 


-0.162 


291.5 


3^Si 


3.200 


0.238 


295.4 


3^ Si 


3.201 


-0.181 


294.8 


3® Si 


3.218 


0.281 


299.8 


3** Si 


3.219 


-0.204 


298.8 


38 Si 


3.224 


0.263 


302.4 


3''Si 


3.245 


-0.254 


301.9 


"°Si 


3.232 


0.244 


305.4 


40 Si 


3.272 


-0.301 


306.0 


"^Si 


3.230 


0.167 


307.1 


"iSi 


3.295 


-0.336 


310.1 


"^Si 


3.228 


0.013 


309.8 


"2 Si 


3.318 


-0.369 


314.6 


"381 


3.240 


0.123 


311.8 


"3 Si 


3.320 


-0.356 


315.2 


""Si 


3.252 


0.172 


314.3 


""Si 


3.322 


-0.342 


316.2 


"5 Si 


3.252 


0.117 


315.8 


"5 Si 


3.316 


-0.308 


317.5 


46 s J 


3.253 


0.053 


317.9 


46 Si 


3.303 


-0.262 


319.3 


"^Si 


3.258 


0.005 


319.7 


"^Si 


3.345 


-0.298 


319.8 


48 Si 


^ 96^ 


n nni 

U.UUl 




48 Si 


J. JO i 


-fl ^9 1 


^90 8 


49si 


3.290 


0.045 


321.1 


49si 


3.366 


-0.251 


320.7 


50Si 


3.319 


0.074 


321.1 


50Si 


3.341 


-0.159 


321.5 


''^Si 


3.345 


0.078 


321.1 


51Si 


3.358 


-0.135 


321.2 


52Si 


3.371 


0.082 


321.4 


52Si 


3.377 


-0.112 


321.2 


«3Si 


3.391 


0.042 


321.6 


53 Si 


3.391 


-0.052 


321.3 


«"Si 


3.415 


0.000 


322.3 


5"Si 


3.415 


-0.010 


322.0 



TABLE V: Same as Table III, but for S 



nucleus 


Vch 


132 


BE(MeV) 


nucleus 


Teh 


B2 


BE(MeV) 


33 s 


3.241 


0.197 


275.5 


33s 


3.233 


-0.116 


275.1 


3"S 


3.248 


0.140 


285.8 


3"S 


3.257 


-0.168 


286.5 


35 s 


3.260 


0.077 


295.6 


35 s 


3.260 


-0.078 


295.7 


36 s 


3.273 


0.002 


306.2 


36s 


3.309 


-0.308 


299.7 


37s 


3.285 


0.152 


311.6 


37s 


3.287 


-0.116 


310.1 


38 s 


3.300 


0.228 


318.6 


38s 


3.300 


-0.164 


316.9 


39 s 


3.312 


0.264 


325.3 


39 s 


3.307 


-0.173 


322.6 


40 s 


3.325 


0.299 


332.4 


40 s 


3.316 


-0.181 


328.5 


41s 


3.331 


0.287 


337.7 


41s 


3.324 


-0.189 


333.6 


42 s 


3.338 


0.277 


343.2 


42 s 


3.335 


-0.207 


339.2 


43 s 


3.359 


0.318 


347.2 


43 s 


3.348 


-0.229 


344.1 


""s 


3.381 


0.367 


351.0 


""s 


3.366 


-0.263 


349.5 


45 s 


3.375 


0.312 


353.4 


45 s 


3.367 


-0.240 


351.8 


46 s 


3.371 


0.258 


356.6 


46s 


3.375 


-0.237 


355.1 


47s 


3.385 


0.257 


358.5 


47s 


3.380 


-0.230 


358.0 


48 s 


3.400 


0.259 


360.8 


48s 


3.389 


-0.230 


360.8 


49 s 




0.227 


^69 Q 


49 s 


3.420 


-0.257 


^69 Q 


50 s 


3.403 


0.189 


365.5 


50s 


3.451 


-0.277 


365.0 


51s 


3.427 


0.188 


366.4 


51s 


3.451 


-0.231 


365.9 


52 s 


3.451 


0.183 


367.6 


52 s 


3.447 


-0.178 


367.6 


53 s 


3.463 


0.158 


369.1 


53s 


3.466 


-0.172 


368.4 


5"S 


3.477 


0.139 


371.0 


5"S 


3.486 


-0.142 


369.8 


55s 


3.494 


0.105 


371.4 


55s 


3.497 


-0.090 


370.5 



this calculations the newly discovered nuclei "Mg and Al 
are well within the prediction both in the SHF and RMF for- 
malisms. Again a further comparison of the drip-line with 
RMF and SHF prediction, we find the drip-line predictions in 
both calculations are well comparable, except for a few ex- 
ceptions in Na and Si as shown in Table 11. 



C. Neutron configuration 

Analysing the neutron configuration for these exotic nuclei, 
we notice that, for lighter isotopes of Ne, Na, Mg, Al, Si 
and S the oscillator shell Nose = 3 is empty. However, the 
Nose = 3 shell gets occupied gradually with increase of neu- 
tron number. In case of Na, Nose = 3 starts filling up at ^^Na 
with quadrupole moment deformation parameter /?2 — 0.356 
and -0.179 with occupied orbits [330]l/2- and [303]7/2-, 
respectively. The filling of Nose = 3 goes on increasing 
for Na with neutron number and it is [330] 1/2", [310] 1/2", 
[321]3/2- and [312]5/2- at P2 = 0-472 for ^^Na. Again for 
the oblate solution the occupation is [301]l/2~, [301]3/2~, 
[303] 5/2- and [303] 7/2" for /32 = -0.375 for ^^Na. In 
the case of Mg isotopes, even for ^°'^^Mg, the Nose = 3 
shell have some occupation for the low-lying excited states 
near the Fermi surface for ^°Mg (at {3= 0.599 with Be = 
237.721 MeV the Nose=3 orbit is [330]7/2- and for ^"^Mg: 
P2 = [330]l/2-, BE=248.804 MeV at ^2 = 0.471). With the 
increase of neutron number in Mg and Si isotopic chain, the 
oscillator shell with Nose = 3 gets occupied more and more. 
For most of the Si isotopes, the oblate solutions are the dom- 
inating ones than the low-lying prolate superdeformed states. 



TABLE VI: The calculated value of charge radii (rch), quadrupole 
deformation parameter (32 and binding energy (BE) for Ne, Mg, Si 
and S even-even nuclei in SHF ( SkI4) formalism. The maximum bind- 
ing energy is the ground state solution and all other values are the 
intrinsic excited state solution. The radius Vch is infm and the bind- 
ing energy is in MeV. 



nucleus 


rch 


Q 

P2 


BE(MeV) 


rch 




P2 


i3ii(Me V) 


JNe 


3.029 


0.5481 


156 


817 


2.950 


-0.1356 


154.474 




3.005 


0.5223 


175 


758 


2.943 


-0.1989 


172.758 


JNe 


2.952 


0.2546 


188 


354 


2.951 


-0.2541 


100 CIO 

188.538 


26-M 

Ne 


2.953 


0.1233 


199 


380 


2.944 


0.0060 


199.389 


Ne 


3.013 


0.1623 


206 


524 


3.010 


-0.1334 


206.433 


3O1VT 
Ne 


3.054 


0.0030 


213 


721 








''^Ne 


3.103 


0.3808 


213 


118 


3.118 


0.3716 


213.215 


Ne 


3.179 


0.4880 


213 


483 


3.108 


-0.1462 


"^/^rv /TIC 

209.695 


Ne 


3.203 


0.6015 


212 


230 


3.147 


-0.2789 


208.770 


Mg 


3.128 


0.5248 


195 


174 


3.077 


-0.252 


ion f\ A /' 

189.946 


26« j„ 
Mg 


3.090 


0.3623 


212 


885 


3.079 


-0.2988 


213.153 


28» j„ 
Mg 


3.111 


0.3419 


228 


997 


3.056 


-0.1076 


227.899 


Mg 


3.119 


0.2022 


240 


328 


3.1 17 


-0.1835 


240.514 


^^Mg 


3.145 


f\ r\f\f\f\ 

0.0000 


252 


033 


3.145 


0.0000 


252.033 


34« T„ 
Mg 


3.209 


0.3263 


255 


067 


3.175 


-0.1 196 


253.455 


^^Mg 


3.295 


0.4884 


201 


512 


3.252 


-0.289 


257.634 


36a t„ 
Mg 


3.265 


0.4413 


259 


899 


3.213 


-0.2124 


255.368 


40« j„ 

Mg 


3.321 


0.4741 


262 


796 


3.299 


-0.3538 


260.200 


28c; 


3. 117 


0.009 


231 


037 


3.194 


-0.3494 


■^o 1 c c\r\ 

233.590 


30e; 


3.145 


0.1477 


252 


146 


3.168 


-0.2102 


252.625 


32c; 
SI 


3.179 


0.007 


269 


479 


3.199 


-0.1990 


270.483 


34 c; 


3.216 


0.000 


286 


332 








36 c; 
SI 


3.146 


0.1549 


292 


418 


3.241 


-0.009 


292.425 


38c; 
61 


3.291 


0.3051 


298 


173 


3.279 


-0.1978 


298.173 


40 c; 
SI 


3.325 


0.2990 


230 


450 


3.309 


-0.2817 


303.969 


42 c; 
SI 


3.349 


0.3592 


307 


399 


3.334 


-0.3508 


310.023 


44 c; 
SI 


3.334 


0.21 19 


309 


712 


3.377 


-0.3031 


311.601 


46 c; 
SI 


3.337 


0.009 


312 


451 


3.372 


-0.2405 


313.508 


48 c; 
SI 


3.348 


0.002 


315 


425 


3.438 


-0.2893 


111 r\c\ c 

313.995 


30 c 
S 


3.262 


0.1491 


241 


434 


3.178 


-0.1856 


241.385 


32 c 
S 


3.271 


0.200 




1 I'X 


3.256 


-0.1700 


267.975 


34s 


3.288 


0.1212 


288 


804 


3.295 


-0.1566 


289.304 


36 g 










3.314 


-0.003 


309.619 


38 s 


3.341 


0.2144 


320 


168 


3.331 


-0.1289 


318.951 


40 s 


3.374 


0.3042 


332 


097 


3.350 


-0.1529 


327.809 


42 s 


3.392 


0.2898 


341 


033 


3.379 


-0.2195 


337.031 


44s 


3.436 


0.3677 


348 


266 


3.410 


-0.2714 


346.445 


46 s 


3.423 


0.2517 


352 


486 


3.413 


-0.2090 


351.578 


48 s 


3.450 


0.2379 


356 


188 


3.430 


-0.2032 


356.589 


50 s 


3.441 


0.1229 


360 


815 


3.498 


-0.2673 


359.011 


52 s 


3.482 


0.1099 


362 


347 


3.487 


-0.1356 


362.531 


54s 


3.528 


0.003 


364 


650 


3.524 


-0.1023 


363.926 


56 s 


3.558 


0.1105 


366 


031 


3.556 


-0.0100 


366.033 



i.e. mass of the oblate solutions are the ground state solutions 
and the prolate and some superdefomed are the excited con- 
figurations. Again, in S-isotopes, the prolate are the ground 
state and the oblate are the extreme excited states. Note that 
in many cases, there exist low laying superdeformed states 
and some of them are listed in the Tables. 




mass number (A) 

FIG. 1 : The ground state quadrupole deformation parameter P2 ver- 
sus mass number A for Ne, Na, Mg, Al, Si and S isotopes near the 
drip-line with NL3 parameter set. 



D. Quadrupole deformation 

The ground and low-lying excited state deformation sys- 
tematics for some of the representative nuclei for Ne, Na, 
Mg, Al, Si and S are analysed. In Fig. 1, the ground state 
quadrupole deformation parameter /32 is shown as a function 
of mass number for Ne, Na, Mg, Al, Si and S. The /32 value 
goes on increasing with mass number for Ne, Na and Mg 
isotopes near the drip-line. The calculated quadrupole defor- 
mation parameter (32 for ■^''Mg is 0.59 which compares well 
with the recent experimental measurement of Iwasaki et al Ht] 
{^2 — 0.58 ± 6). Note that this superdeformed states in 3.2 
MeV above than the ground band. Again, the magnitude of 
/32 for the drip-line nuclei reduces with neutron number N and 
again increases. A region of maximum deformation is found 
for almost all the nuclei as shown in the figure. It so happens 
in cases like, Ne, Na, Mg and Al that the isotopes are maxi- 
mum deformed which may be comparabled to superdeformed 
near the drip-line. For Si isotopes, in general, we find oblate 
solution in the ground configurations. In many of the cases, 
the low-lying superdeformed configuration are clearly visible 
and some of them are available in the Tables. 



E. Shape coexistence 

One of the most interesting phenomena in nuclear structure 
physics is the shape coexistence ll29l - [3Tll . In many of the cases 
for the nuclei considered here near the drip-line isotopes, the 
ground state configuration accompanies a low-lying excited 
state. In few cases, it so happens that these two solutions 
are almost degenerate. That means we predict almost similar 
binding energy for two different configurations. For example, 
in the RMF calculation, the ground state binding energy of 
24Ne is 189.093 MeV with P2 = -0.259 and the binding en- 
ergy of the excited low-lying configuration at P2 = 0.278 is 
188.914 MeV. The difference in BE of these two solutions is 
only 0. 179 MeV. Similarly the solution of prolate-oblate bind- 
ing energy difference in SkI4 is 0.186 MeV for ■^°Mg with 



■ I ■ ■ I ■ ■ I ■ ■ I ■ ■ I ■ ■ I 



I ■ ■ I ■ ■ I 



24 30 36 42 48 54 6024 30 36 42 48 54 60 
A 

FIG. 2: The difference in binding energy between the prolate-oblate 
solutions is shown for even-even Ne. Me. Si and S isotoues near the 



50 — neutron 



P=0.S88 



proton 



proton 
|3=0.343 



proton proton 
M-471 M-119 



— I — ' — r 

SHF (SkI4) 




20 24 28 32 36 4020 24 28 32 36 40 
N 

FIG. 3: The two-neutron separation S2n energy versus neutron num- 
ber N for neutron-rich Ne, Mg, Si and S isotopes. 



^2 = —0.183 and 0.202. This phenomenon is clearly avail- 
able in most of the isotopes near the drip-line. To show it in 
a more quantitative way, we have plotted the prolate-oblate 
binding energy difference in Figure 2. The left hand side of 
the figure is for relativistic and the right side is the nonrela- 
tivistic SkI4 results. From the figure, it is clear that an island 
of shape coexistence isotopes are available for Mg and Si iso- 
topes. These shape coexistence solutions are predicted taking 
into account the intrinsic binding energy. However the actual 
quantitative energy difference of ground and excited config- 
uration can be given by performing the angular momentum 
projection, which is be an interesting problem for future. 



F. Two neutron separation energy and new magic number 

The appearance of new and the disappearance of known 
magic number near the neutron drip-line is a well discussed 
topic currently in nuclear structure physics fs', 'sS]. Some of 
the calculations in recent past predicted the disappearance of 
the known mag ic number N=28 for the drip-line isotopes of 
Mg and S l33l l34ll . However, magic number 20 retains its 



FIG. 4: The l/2^andl/2^ intrinsic single-particle states for the 
normal and superdeformed state for ^^Mg and '^'*Mg. A few of the 
lowest energy parity-doublet states of the superdeformed (SD) so- 
lutions are shown by asterisk for the SD configuration. More such 
doublets are noticed for the SD intrinsic states. The ±l/2~ states 
are denoted by shorter (and green) lines and the ±1/2"'" states are 
denoted by longer ( and black). 



magic properties even for the drip-line region. In one of our 
earlier publications, ifssll we analysed the spherical shell gap 
at N=28 in ^''S and its neighboring ^°Mg and ''^Si using NL- 
SH ||22ll and TM2 parameter sets [36]. The spherical shell 
gap at N=28 in ^*S was found to be intact for the TM2 and 
is broken for NL-SH parametrization. Here, we plot the two- 
neutron separation energy 82,1 of Ne, Mg, Si and S for the 
even-even nuclei near the drip-line (fig 3). The known magic 
number N=28 is noticed to be absent in *^S. On the other hand 
the appearance of steep 2n-separation energy at N=34 both in 
RMF and SHF calculation looks quite prominent, and this is 
just two units ahead than the experimental shell closure N=32 



G. Superdeformation and Low Q parity doublets 

The deformation-driving m = 1/2— orbits come down in 
energy in superdeformed solutions from the shell above, in 
contrast to the normal deformed solutions. The occurrence of 
approximate 1/2+ 1/2^ parity dobulets (degeneracy of |m|'^= 
1/2+ 1/2^ states) for the superdeformed solutions are clearly 
seen in Figs. 4 and 5 where excited superdeformed configura- 
tions for •^•^Mg and ■^''Mg and for '*^A1 and ''^Al are given. 
For each nucleus we have compared the normal deformed 
(/?2 ^ 0.1 — 0.3) and the superdeformed configurations and 
analysed the deformed orbits. 

The 1/2+ and l/2~ states for the single particle levels are 
shown in Fig. 4 (for '^^Mg and "^^Mg). From the analysis 
of the results of this calculation, we have found a systematic 
behaviour of the low fl (particularly 1/2+ and 1/2^) prolate 
deformed orbits for the superdeformed solutions. As repre- 
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Al (1/2^ & 1/2") Al 







~ neutron 

proton proton 
- M-653 [3-0.090 


neutron 

neutron — 

M.660 " 



FIG. 5: Same as Fig. 4 for -^^Al and Al. 



sentive cases, we present here results for Mg — ^^Mg) and 
(47^1 _ 46^1) jjjj^j piQj jijg 112+ and 1/2- orbits for the su- 

perdeformed and normal deformed shapes of these nuclei. We 
notice from the plot of the orbits that there is occurrence of 
1/2+ and 1/2^ orbits very closeby in energy for the superde- 
formed (SD) shape. Two such 1/2+ and 1/2^ doublet struc- 
tures, marked in asterisk are shown in Figs. 4 and 5 for the 
SD solutions. Such 1/2+ 1/2^ degenerate orbits occur not 
only for the well-bound orbits but also for the unbound con- 
tinuum states. As example, the doublet neutrons [220] 1/2+ 
and [101]l/2^ states is 4 MeV apart in energy in the normal 
deformed prolate solutions tend to become degenerate in the 
SD solution [220] 1/2+ and [101] 1/2" states (prolate) belong- 
ing to two different major shells, so close to each other in the 
superdeformed solution (shown in Figs. 4 for "^"^Mg). More 
such doublets are easily identified (Figs. 4 and 5) for superde- 
formed solutions of ^^•■^''Mg and ^^■^''Al. In fact it is noticed 
that the = 1/2 states of unique parity, seen clearly well 
separated in the normal deformed solutions, get quite close to 
each other for the SD states, suggesting degenerate parity dou- 
blet structure. This can lead to parity mixing and octupole de- 
formed shapes for the SD structures jssll . Parity doublets and 
octupole deformation for superdeformed solutions have been 



discussed for neutron-rich Ba and Zr nuclei There 
is much interest for the experimental study of the spectra of 
neutron-rich nuclei in this mass region 140]. The highly de- 
formed structures for the neutron-rich iVe — Na — Mg — Al 
nuclei are interesting and signature of such superdeformed 
configurations should be looked for. 



IV. SUMMARY AND CONCLUSION 

In summary, we calculate the ground and low-lying excited 
state properties, like binding energy and quadrupole deforma- 
tion /32 using NL3 parameter set for Ne, Na, Mg, Si and S 
isotopes, near the neutron drip-line region. In general, we find 
large deformed solutions for the neutron-drip nuclei which 
agree well with the experimental measurement. We have done 
the calculation using the nonrelativistic Hartree-Fock formal- 
ism with Skyrme interaction SkI4. Both the relativistic and 
non-relativistic results were found comparable to each other 
for the considered mass region. In the present calculations 
a large number of low-lying intrinsic superdeformed excited 
states are observed for many of the isotopes and some of them 
are reported. The breaking of N=28 magic number and the 
appearance of a new magic number at N=34 appears in our 
calculations. A proper angular momentum projection may tell 
us the exact lowering of binding energy and it may happen 
that the superdeformed would be the ground band of some of 
the neutron-rich nuclei. Work in this direction is worth doing 
because of the present interest in the topic of the drip-line nu- 
clei. In this study we find that, for the SD shape, the low VI 
orbits (particularly = 1/2) become more bound and show a 
parity doublet structure. Closelying parity-doublet band struc- 
tures and enhanced electromagnetic transition rates are a clear 
possibility for the superdeformed shapes. 
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